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Abstract 

We construct a definition of the weak solution to KdV type equations with small 
dispersion admitting the zero dispersion limit for soliton-like solutions. Using this 
definition, we obtain a system of equations (the limit problem as the dispersion tends 
to zero) that describes the soliton dynamics. 



1 Introduction and basic results 

1 . It is well known that the Korteweg-de Vries (KdV) equation 

L K dv[u] = u t + {u 2 ) x + s 2 u xxx = (1.1) 

has the one-soliton solution 

u(x,t,e) = ^ch- 2 (^(x-vt)/e), ieR, (1.2) 



where v is the soliton velocity. The pointwise limit as e — > +0 of solution ( |l.2| ) to the KdV equation 
is the discon tinu ous function ^x( x — where %(£) = 1, if £ = and y(£) = if £ 7^ 0. The weak 
asymptotics ( |l.2| ) as e — > +0, up to Op'(£ 2 ), becomes the infinitely narrow 5-soliton 

u E (x,t) = AeS{x - vt), A = — I ch~ 2 (£) d£ = 6y/v, (1.3) 

and <5(x) is the Dirac delta function. Here and in what follows J denotes an improper integral from —00 
to +00. By 0-u>{e a ) we denote a distribution from Z>'(R) such that for any test function tp(x) G V 

(O v >(e a ) > <p(x))=0(s a ), 

and 0(e a ) is understood in the ordinary sense. 

We stress once more that here all generalized functions (distributions) are treated as functionals on 
the space T>(EL X ) and these functionals depend on the other variables as on parameters. 



It follows from (1.2), (1.3) that we have u(x,t,s) — Ot>'(e) as e — > +0 and e 2 u xxx = Ox>/(e 3 ). There- 
fore, the limit expression (1.3) was interpreted by V. P. Maslov and V. A. Tsupin V. P. Maslov 
and G. A. Omelyanov g-gf as an asymptotic up to Ot>>(s 2 ) generalized solution of the Hopf equation 

L H [u] =u t + {u 2 ) x =0, (1.4) 

which is the limit problem for the KdV equation. In the same works the corresponding generalized 
Hugoniot conditions, of the type of those for the shock wave front, were obtained. 
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It is easy to see that in the framework of the above-mentioned approach, the exact solution (1.2) to 
KdV equation (IT) satisfies the Hopf equation (L4) in the following sense: 



L H [u] = O v >{e 2 ), 



because the following equality holds: 



LKdv[u] - L H [u] = Ov{e 2 ). 



(1.5) 



(1.6) 



If, instead of expression (1.2) which is the exact solution of the KdV equation and approximates the 
weak asymptotics (|l.3D, we consider the function 



u(x,t,e) — Auj( ), 

e 

has a compact suppor t or rapidly decreases as \z\ 



(1.7) 

* oo, / u){z) dz = 1, then: 
+0 (see Section || for details); 



where uo(z) £ C° 

1) expression (1.7) also has asymptotics (L3) i n th e sense of V as e 

2) substituting (1.7) into the Hopf equation (1.4), provided that a certain correlation between the 
constants v and A is true (the generalized Hugoniot condition), we have = Ot>>(s 2 ); 

3) and u(x, t, e) — u(x, t, e) = Ox» (e 2 )- 

Therefore, an asymptotics up to Ot>'(s 2 ), i.e., an infinitely narrow (5-soliton-type solution of the KdV 
equation (or the Hopf equation which is the limit problem of the KdV equation) , can be sought starting 
not from the exact solution (1.2) of the KdV equation (which is a regularization of the Hopf equation) 
but from an ansatz of the form (1.7) substituted directly into the Hopf equation. 

Generalizing (1.3), one can seek the solution in the form 



4>(ty 



u*(x,t,e) = u (x,t) + g(t)u(- 

by subsituting this singular ansatz into the Hopf equation. 

However, as we shall see in Section [|l, such an ansatz results in the solution with constant amplitude 
g = const of the soliton for u$ ^ const, which contradicts the well-known results about the soliton 
behaviour. Therefore, generalizing formulae ([^), ( |1.7| ), ( |1.3j ), we can attempt to construct an asymptotic 
solution to the KdV equation (or the Hopf equation) of the form 



u*{x, t) — u (x, t) + g{t)e5(x — (f>(t)) + e(x, t)eO(x — (f>(t)), 



+0, 



(1.8) 



where 6 is the Heaviside function. A solution in such a form will be called an infinitely narrow S-soliton. 

In order to consider the function u*(x,t) as a solution to a nonlinear equation we need to define the 
rules of substitution of this function into the nonlinear term, that is, to include the functions e5(x — <j)(t), 
e9(x — 4>(t)), where e — > +0, into an algebra with differentiation. 

2. The rules of substitution of singular ansatzs into nonlinear equations are finally reduced to the 
definition of multiplication of distributions from certain classes. Various approaches to this problem are 
discussed, for example, in [^J . A well-known approach to the problem of multiplication of distributions 
related to the name of J. Colombeau 0, M is based on the construction of an algebra of new generalized 
functions. The spaces C°° and D 1 are embedded into G, and moreover, the first embedding is an algebraic 
isomorphism. Some of the new generalized functions are associated with distributions but, in the general 
case, they are not directly related to the Schwartz distributions. 

Note that H. A. Biagioni and M. Oberguggenberger (S) studied the solution to the KdV equation 
in the form of an infinitely narrow soliton (1.2) in terms of the theory of Colombeau new generalized 
functions. 

Another approach to this problem is associated with the names of J. B. Livchak §, Li Bang-He 0, 
V. K. Ivanov ]H[ . The main idea of this approach is that products of distributions from a certain class 
are defined as asymptotic decompositions, in the weak sense, of products of approximations of these 
distributions as e — > +0, where e is an approximation parameter, and the coefficients are distributions. 

Development of ideas of this approach in , Jlj| and especially in [jl4| has led to the construction of 
the weak asymptotics method which enables one to obtain substantial analytical results for discontinuous 
solutions to nonlinear equations and, in particular, to study the dynamics of propagation and interaction 
of singularities of various types. In the framework of this approach the construction of a singular solution 
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to a nonlinear equation is reduced to the choice of singular generators of the ansatz, that is, distributions 
(generalized functions), to defining the rules of substitution the solution with the chosen structure into the 
nonlinear equation, to setting up a system of equations for unknown smooth functions and to investigation 
of the system obtained. 

An essential technical progress achieved in these papers, as compared with the results of |J] - |rH| , is 
the obtaining of complete asymptotic decompositions of products of approximations of distributions and, 
on this basis, of associative and commutative differential algebras of asymptotic distributions including 
subspaces of distributions. The elements of these algebras are weak asymptotics whose coefficients are 
distributions. Thus, in Jl4| there was constructed an associative and commutative algebra £ of asymptotic 
distributions (weak asymptotics) generated by the space of linear combinations of homogeneous or asso- 
ciated homogeneous distributions. For example, the singular ansatz (1.8) is an asymptotic distribution 
belonging to this algebra. 



An asymptotic generalized infinitely narrow <5-soliton-type solution to the Hopf equation (1.2) with 
constant amplitude, and the generalized Hugoniot condition were constructed in the framework of the 
technique of substitution of a singular ansatz into [|12l , flU 




The application of the weak asymptotics method in ]14| has made it possible to obtain some substantial 
results about the structure of singular solutions to quasilinear strictly hyperbolic systems and, with 
certain restrictions, to find the solution of the problem of separating self-similar singularities formulated 

in una, M. 



3. Th us, in order to substitute the singular ansatz (L8) into the KdV equation (IT) or the Hopf 
equation (IT), we first substitute there a smooth ansatz of the form 



u*(x, t, e) — u (x, t) + g(t)eS(x — <f){t), e) + e(x, t)e6(x — (j)(t), e), e > 0, 
where Uq(x, t), g(t), e(x,t), <j){t) are the desired smooth functions, and 

e6(x, e) = u(— ), e9(x,e) — eujo( — ) 



(1.9) 



are smooth approximations of the asymptotic distributions e5(x) and s9(x), respectively. 

Here, the function oj(z) S C°°(R) either has a compact support or decreases sufficiently rapidly as 
\z\ — > oo, for example, \ui(z)\ < C(l + \z\)~ 3 and J uj(z)dz = 1; ui Q (z) E C°°(R), lim z ^ +oo ui (z) = 1, 
lim^-oo lo (z) = 0. 

Then we have in the sense of 2?'(IR) (see the notation above) 



e5(x,e) = eS(x) + O v ,(e 2 ), e6(x,e) = e9{x) + B »(e 2 ), 



-0. 



For more details, see Section S 



Now, by analogy to (1.5), (1.6), we can introduce the definition of the asymptotic generalized solution 
of the form (1.?). Namely, we call asymptotic distribution fll.SP an asymptotic generalized solution of the 



Hopf equation (1.9) if its approximation (l.£) satisfies the relation 

L H [u*(x,t,e)] = O v >{e 2 ) 



(1.10) 



and, respectively, an asymptotic generalized solution to the KdV equation (1.1), if 

L Kdv [u*(x,t,e)] = Ovi(e 2 ), (1.11) 

which, in fact, is the same, since 

L KdV [u*(x,t,e)} - Lff[u*(x,t,e)] = O v >{e 2 ). 

It is easiy to see that our definition of the solution can depend on the choice of approximations 
-a>( z ~'?ty ) and ujo(— to the distributions S(x — 4>(t)), and 9(x — <ft(t)), respectively. Actually, as 
we shall see later, the dynamics of solution of the type (1.8) is independent of the approximation of the 
Heaviside funct ion u) o( x ~f ^ -). 

Definitions ( l-lCQ and jLll|) als o imply that, in fact, in order to construct asymptotic generalized 
solutions satisfying ( l.lCj ) or (1.11), we need not calculate products of generalized functions (and, in 
general, any other nonlinearities) to a high accuracy. Indeed, the substitution of the exact one-soliton 
solution of the KdV equation into the Hopf or KdV equations provides just the accuracy (0-p/(e 2 )) 
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corresponding to ( 1.10 ). On the other hand, it is clear that we cannot deal with a lesser accuracy, since 
the soliton solution contains terms of order Ox>>(s). 

In our notation, we can say that, within the framework of algebraic constructions related to general- 
ized functions, nonlinear expressions are usually calculated up to Od»(£°°), which is necessary to define 
associative and commutative algebras of generalized functions. 

A distinction of the method of weak asymptotics (which, undoubtedly, originates from the algebraic 
constructions mentioned above) is that we actually deal with approximations. In fact, the difference 
between the method of weak asymptotics and the method of ordinary asymptotic expansions is that the 
smallness of the remainder is understood in a different way. Usually, the remainder is assumed to be 
small in some uniform sense with sufficient accuracy. Here we assume exactly the same but in the sense 
of Ot>'. 

In order to obtain the results known from the KdV equation theory it seems natural to use the function 



from the formula for the exact one-soliton solution (1.2) to the KdV equation as an approximation for 
e5(x — 0(i), e). 



The system for the functions uo(x,t), g{t), e{x,t), <f>{t) follows from Definitions (1.10) or (1.11) (this 
system will be derived in detail in Section 1.3) 



e(J>(t),t)-^-g t (t)/g 3 / 2 (t) 



(e t (x, t) + 2(u (x, t)e(x, t)) x ) 



x><j>(t) 



0. 
0. 

0. 
0. 



(1.12) 



It is easy to verify that under the condition g > (which is an analog of the admissibility condition 
in the theory of shock waves) the solution of system ( 1.1 2| ) exists on any interval t G [0, T] such that the 
smooth solution uq of the Hopf equations exists on this interval. 

System (|TT|) can be solved in the following way: first, one finds the smooth solution of the Hopf 
equation, next, one finds the function e(x, t) from the last equations (which is uniquely solvable in view 
of the inequality 2uo(cf),t) < <f> t ), then one finds the (positive) function g(t) from the next to the last 
equation, and finally, one finds the function <f>{t). 

Note that system (1.12) contains no obstacles to setting e(x,t) = 0. If so, g(t) — const in the case of 
an arbitrary (nonconstant) background function uo(x,t). But this conclusion is contrary to well known 
properties of soliton solutions of the KdV equation (see, e.g., ||). 

Moreover, under our notation, the weak asymptotics of the asymptotic one-soliton solution to the 
KdV equation, constructed by V. P. Maslov and G. A. Omelyanov has the form 



u{ £ {x, t) = u i(x, t) + gi {t)e5{x - 0i (t)) + ei (x, t)e[l ~ 6(x - <M*))]> 



-0. 



(1.13) 



In other words, in the case (1.8) the "shock wave" with a small amplitude se(x, t)6{x — 4>i(t)) propagates 
in front of the soliton eS(x — (f>i(t)), but in the asymptotic one-soliton solution constructed in |j| the 
small shock wave eei(x,t)[l — 9{x — <fii(t))] arises behind the soliton . 

If we apply Definition (1.10) or (1.11) to the asymptotic solution obtained in ||, whose weak asymp- 
totics yields ( 1.13Q , we obtain the following system of equations 



"Olt + («oi)a: 

0it-2«oi(0 1 (t),t)- - gi {t) 
eim),t) + ^9it{t)/gl'\t) 



(e u (x,t) + 2(u i(x,t)e 1 (x,t)) x ) 



x<4>i(t) 



0. 
0. 

0. 
0. 



(1.14) 



The solution of the last system for guit) ^= is not uniquely determined by the initial conditions 
ei(x, 0) for x < (j>i(0), since the velocity along the characteristic (x = 2uoi(x(t),t)) is less (for gi(t) > 0) 
than that of the soliton 4>it = 2u i((j)i(t), t) + §<7i(t). 



Thus, the assumption that the structure of the solution to the KdV equation is specified by ( |1.13| ) 
due to Definitions fll-lCQ , (1.11) leads to an ill-posed Cauchy problem (with a nonunique solution) for the 
functions Mqi (x, t), gi(t), ei(x,t), 4>\{t). 
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On the other hand, the system of equations obtained in B for these functions has the form 



uqu 4 




= o, 


<f>it - 2«oi (0i (*)>*) - 




= o, 


ei(<K*),*) + ^W*)/s? /2 (*) 


= o, 


+ 2(uoi(x,t)ex(x,t)) x ) 
gi(t) + 2u i((, 


x<4>i(t) 
il(t),t) 


= o, 

= const, 



(1.15) 



It is evident that this system differs from system (1.14) obtained from (1.1C), (1.11), (1-13) by the 
additional eq uatio n gi(t) + 2uoi(4>i(t), t) = gi(0) + 2uqi(0i (0), 0). The presence of this equation implies 
that system ( 1.15 ) splits into the two systems 



2uoi(^i(t),«)-^i(*) 
9l (t) + 2u 01 (Mt),t) 



0, 

0, 

const, 



and 



(e lt (x,t) + 2(u 01 (x,t)e 1 (x,t)) x 
3\/6 



eiW),t) + ^ r g lt (t)/gt''{t) = 



x<<t>i(t) 

3 / 2 n\ - 



0. 



(1.16) 

(1.17) 
(1.18) 



and equality ( 1.1 8| ) is the boundary condition for equation (1.17), which turns the Cauchy problem for 
equation ( 1.17| ) into the well-posed one (the Cauchy condition, in view of ( 1.13| ), has the form e%(x, 0) = 

4(x)[i-0(x-moM 

Moreover, if equation ( |l,17 ) is considered formally in the domain x > <f>x (t) , which corresponds to the 
solution structure given by formula (1.8), then the "redundant" condition 



e(4>(t),t) 



3^6 



9t(t)/g 3/2 (t) = 0, 



analogous to (1.18), overdetermines the problem. 

Thus, the weak asymptotics corresponding to the asymptotic solution of the Cauchy problem for the 
KdV equation const ruct ed in [|| cannot be derived from the solution to the KdV equation with the help 
of Definitions ( 1.10 ) or ( 1.11 ), and vice versa. 

Why is it so? The essence of the matter lies in the definition of weak (generalized) solution to nonlinear 
equation. It turns out that the definition of the weak (generalized) solution to nonlinear equation depends 
on the structure of the kernel of the operator adjoint to the linearized operator of the initial differential 
equation which arises when constructing the smooth asymptotics. This construction of the definition of 
weak solutions was previously discussed in Q , Q . 

In the present paper we do not come into details of construction of the definition of the weak solution 
to our problem. We just point out that, in terms of this construction, the KdV equation is analogous to 
the phase field system discussed in 

The difference is that for the KdV equation the kernel of the adjoint operator mentioned above is 
two-dimensional, which results in the following definition analogous to that of the weak solution in the 
form of the integral identity from |Q . 

Definition 1.1 The asymptotic distribution of the form ( 1.13| ) 

u* £ {x,t) = u (x,t) + g{t)e8(x - <f>(t)) + e(x,t)s0(-x + <f>(t)) 

is a generalized asymptotic (soliton-type) solution to the KdV equation ([0]) for t S [0, T] with the initial 
condition u®*(x), if for any constants c%, C2 the following equality holds 

(ci +c 2 u*(x > t,e))L Kd v[u*(x,t,e)] = O v >{e 2 ), 

<(x,0) = u E *(x)+O T , l (e 2 ), 



(1.19) 



where u*(x, t, e) is a smooth approximation of the asymptotic distribution u*(x, t), and the first estimate 
is uniform with respect to t € [0, T]. 
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It is clear that ( 1.19 ) is equivalent to the following relations 

LKdv[u*{x,t,e)) =Ov(e 2 ), u*{x,t,e)L KdV [u*(x,t,e)} = O v >(e 2 ), 



and the first relation coinci des w ith Definition ( 1.10 ). 

One can easily see that ( 1-1S| ) can be rewritten as an integral identity but of an unusual form. 

By analogy to what was previously said, the solution depends on the choice of the approximation, and 
to obtain the results known in the theory of the KdV equation one should choose, as an approximation 
of the asymptotic distribution g(t)eS(x — </>(i)), the function from the formula for asymptotic solution to 
the KdV equation ||: 

g(t)eS(x - cf>(t),e) = g(t)u(a(t f ~ ^ ' ), 

where a(t) = *J ^[p, oj{z) = ch~ 2 (z). 

It is well known that the function u>{a{t)r) is a solution of the boundary value problem for the 
differential equation 

, -.dbJ , , , . . . . . . dw o , > <Plo 

-&{t)— + 2{u {cj ) (t),t)+g{t)Lo)— + a 2 {t)—=Q, (1.20) 

where u)(z) — > as \z\ — > oo. 

It is easy to see that this boundary value problem has the solution given above provided that 



<Pt = 2u (<t>(t),t) + ~g(t). 

In fact, the equation for the approximation u> can be derived by applying an analog of the average 
procedure, but in the present paper we shall simply postulate the choice of approximation (1.20). 

In Section ^ it is shown that for approximation of a soliton of the form (1.20), Definition 1.1 imme- 
diately leads to the system of equations (1.15) previously obtained in 
asymptotics of one-soliton solution to the KdV equation. 
It should be emphasized that Definition 



when constructing the smooth 



1.1 



for the Hopf equation is not equivalent to that for the 
KdV equation since the term e 2 u xxx , for C2 ^ 0, makes a contribution of the order of O-pi(e) to the 
left-hand side of the equality. 

Thus, it can be said that asymptotic generalized solutions exist for both the Hopf and KdV equations . 
When ua(x,t) — const, the systems describing this solution coincide but, in the general case (uq(x, t) ^ 
const), they prove to be absolutely distinct in the sense that the limit problem for any one of them is 
ill-posed for the other. 

Moreover, the formal <5-soliton generalized asymptotic solution of the Hopf equation is not related (is 
not an asymptotics) to any exact solution of the Hopf equation. In order to explain this fact, we construct 
the exact solution of this equation satisfying the Cauchy condition 

W(x,t,e)\ t=0 = e5(x). 

For t > this (discontinuous) solution has the form 

0, x<0, 
W(x, t, e) = { x/2t, < x < 
0, x > 2y/ei. 



The discontinuity points </)(t) = 2y'et can be calculated from the Rankinc-Hugoniot conditions. Let 
tp(x) E T>(M. X ) be an arbitrary test function. Let us calculate the limit 



lim {W(x,t, e),ip(x)) = lim 



1 



xip(x) dx. 



By using the L'Hospital rule, we obtain 



\im o (W(x,t,e),<p(x)) =e<p{0). 
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It is clear that for t = 1 the support of W(x,t,e) is the interval [0, 2^/e\, wh ile fo r t = t Q the singu- 
larity support of the generalized (5-soliton constructed, in view of Definition ( jl.lOD , with an arbitrary 
approximation lj(x/e) of the initial value sS can be found at the point 

0(*o) = J u 2 (z) dz, 

which lies outside the support of the exact solution W(x, t, e) (the method used for deriving this relation 
will be discussed below). 

The results, similar to mentioned above can be obtained for the KdV type equation u t + (f(u)) x + 
£ 2 u X xx = with arbitrary smooth nonlinearity, which admits soliton solutions. To verify this assertion, 



by analogy to (1.20), we consider the boundary value for the ordinary differential equation 

, duj d ,, , , . d 3 

-fo-T + l-/( u o(0 ; t)+ui) + —w = 0, 
dr dr dr a 

where uq, 4> are treated as parameters and the solution u> = u>(t, t) is chosen as an approximation of the 
term eg(t)5(x — 4>(tjj in the singular ansatz. Integrating this equation, multiplying the obtained equation 
by oj t , and integrating again, we obtain the energy conservation law 

1 ~(u T f + <t>{u) = E, 

where the function t) = f Q f(z + uq) dz — ^ <?W 2 is considered as a potential and E is an arbitrary 
constant having the meaning of energy. 
Integrating the last equation, we find 

T ~ T Q = ± [" - ^ (1.21) 

where cji, L02 are the last two roots of the equation $(£) = E, the radicand is assumed to be positive 
when u>% < u>(t) < u>2, tq is an arbitrary constant. 

Using the standart met hods, one can show that, for a certain choice of the function f(u) and the 
constant E, equation ( 1.21 ) has soliton solutions (see jl{], §1]). 



For the existence of localized solution it is necessary that one of the roots u>x, u>2 of the equation 



$(£) = E is multiple (that is, the integral in ( 1.21 ) is divergent). 

Thus, in the case of equations with general nonlinearity, we shall substitute into them the smooth 
ansatzs 

x — d>(t) X — d>(t) 

u*{x,t,e) = uo(x,t) +u>( — ,*) +e(x,t)euj ( ^), e > 0, (1.22) 

£ £ 

where the function uj(t, t) decreases sufficienly rapidly with respect to r and luq^x/e) is an approximation 
of the Heaviside function 9(x). 



Accordingly, in this case the weak asymptotics (1.22), that is, for equations of Hopf or KdV type, the 
singular ansatz has the form 

u* e (x,t) = u (x,t) + fl 1 (t)eS(x- <f>(t)) +e(x,t)e0(± (x-(j>{t))), e -> +0, (1.23) 



til(t) = J u{T,t)dT. 

Note that, generally speaking, we need not use any fixed approximations for the distributions contained 



in (1.23). In this case, we obtain equations that describe the dynamics and, in general, depend on the 
approximation of singular generators of the ansatz. 

The dependence of the asymptotic discontinuous solution to a nonlinear equation on approximations 
of singular generators of the solution is substantial for our approach. It is due to the fact that asymptotics 
of the product of approximations is independent of the choice of these approximations only in specific 
cases, for example, when the factors are "weakly" singular, as when calculating the asymptotics of the 
product of the aproximation of a distribution by the approximation of a smooth function. Naturally, 
there is no such dependence if we use this method to find discontinuous solutions of linear equations. 
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By what has been said, when constructing the set of asymptotic distributions in which singular ansatzs 
of the form (1.8), ( 1.13 ) are embedded all possible approximations will be considered for each distribution. 



The same idea is basic to the construction of the Colombeau algebra of new generalized functions || , (?]] 
In the Colombeau theoty there are also "many" (^-functions associated with "one" Schwartz's 6- function. 
Let us summarize the main results. One of the most important results is Definition 1.1 for 



4. 



generalized asymptotic solutions of KdV type equations. In order to explain how to use this defintion, in 
Section || of the present paper we construct an associative algebra of smooth functions tt* in which the 
smooth ansatzs (1.10), (1.14) or (1.22) are embedded which are approximations of the singular ansatzs 
( p^ ), ( fLlEj) or (|L23f ). The singular ansatzs on which generalized infinitely narrow soliton-type solutions 
are constructed belong to the set of asymptotic distributions £* and are derived as weak asymptotics of 
elements from tt* . In this section we also give a definition of generalized solutions to nonlinear equations, 
that is, rules of substituting singular ansatzs into these equations. 

In Section [}] we derive the above-mentioned systems of equations (1.12) and (1.15) by using different 
definitions of the generalized solution. 

In Section |] we write out the system of equations determining the dynamics of the infinitely narrow 
(5-soliton to equations of Hopf and KdV types with general nonlinear terms. 



2 Asymptotic distributions and generalized solutions 
to nonlinear equations 

1. When solving problems listed in Introduction it is necessary to construct a singular solution with 



singular generators of the special type (1.8) or (1.13) to the nonlinear equation L[u] = 0. According to 



our method of weak asymptotics, we substitute into the equation L[u] = smooth ansatzs of the form 



(1.9) resulting from the replacement of singular generators by their approximations in singular ansatzs 
of the form Q. 

It is well known |Hl ch.I, §4.6] that to each distribution (generalized function) f(x) € V can be 
assigned its approximation 

f(x,e) = f(x)*K(x,e) = (f(t),K{x-t,e)), e > 0, (2.24) 

where * is a convolution, the kernel K(x, e) — ^w(f) is a S-type function such that ui(z) € C°°(R), u>{z) 
has a compact support or decreases sufficiently rapidly as \z\ — > oo, for example, \u){z)\ < C(l + |z|)~ 3 
and J uj{z) dz = 1. 

For all test functions <p(x) £ T> we have: 

lim (f(x,e),ip(x)) = (f(x), <p(x)). 

e— »+0 

We cite the approximations of the distributions 5(x), 0(x), which are singular generators of the 
soliton-type ansatzs (1.8), ( 1.13| ). 



For the approximation of 5- function we have from ( |2.24 ) 



8{x,e) = -uj{-). (2.25) 

e e 

For the approximation of the Heaviside function 9(x) we have 

i r°° 
9(x, e) = 6{x) * = / w(- - t) dt. 

Hence we find 

0(a:,e)=wo(~), (2.26) 

where oj (z) = /^^(r/) drj, Yua z ^ +00 u (z) = 1, lim ^_o O ^ (z) = 0, ujq(z) e (7°°(R). 

If the Cauchy kernel K(x,e) — ^ xi e + a is used in ( 2.24 ), we obtain harmonic approximations f(x,e) 
for distributions f(x). In this case (see Q, p2|): 

6( m -V(x,e) = (" 1 )>" 1 ) ! (,-m_ r m )| 6{x ^ £) = fi + l a rctg(-)V 
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where z = x + ie, z = x — ie, m = 0, 1, 2, . . .. 

Denote by TL* the associative and commutative differential algebra generated by finite sums of finite 
products of functions f(x,e) from the space TLq of smooth ansatzs u*(x, t, e) of the form (!.£), where 
approximations of the <5-function and the Heaviside function are defined in (2.25) and (2.26), respectively. 

Denote by £* the set of weak asymptotics f*(x,t) derived from elements f*(x,t,e) of the algebra 
TL*, as e — > +0, determined up to 0-p'(e 2 ), and, according to jyj, [li] , call it the set of asymptotic 
(soliton-type) distributions. 

By analogy with |l4|j , one can introduce a structure of the asymptotic algebra on £* which is not, 
however, required in the present paper. 



Thus, for example, the one-soliton solution to the KdV equation (1.2) belongs to the algebra TL* and 
its weak asymptotics ( |l.3|) a nd singular ansatzs of the form (1.18) are asymptotic distributions from £*. 

It can be shown (see pl4|) that each element f*(x,t,e) from the algebra TL* has a weak asymptotics 
of the form 



f*(x, t, e) = f Q (x, t) + fi(x, t)eS((x - ^t))- 1 ) + f 2 (x, t)e6((x - ^.(t))" 1 ) + O v > {e 2 ), 



-0, (2.27) 



where fj(x, t) are smooth functions, j — 0, 1, 2. 

Finding weak asymptotics of the elements from the algebra TL* can be interpreted, in a sense, as 
multiplications of asymptotic distributions (distributions) (see [ fTpf , Q and |Q- Jl4[). 

Let us give some examples. The following equalities can be readily verified: 



[s6(x,e)] n 
e6(x, e) 
[e9(x,s)] m 
[e5{x 1 e)\ n [ee{x 1 e)] m - 1 



= a n eS(x) + V '(e 2 ), 

= e9(x) + O v >(e 2 ), 

= O v ,(s m ), 

= O v> (e m ). 



(2.28) 



where n = 1, 2, . . ., m — 2, . . ., lu(z) is a <5-type function from ( ^.24 ), fl n — u) n (z) dz. 

In particular, it follows herefrom that the asymptotic distributions sS(x) and ed(x) constitute an 
algebra, which is asymptotic modOx><(e 2 ). 

2. Let F(u) be a smooth function of at most exponential growth. We define a function of asymptotic 
distribution F(f*(x,t)) £ £* as the weak asymptotics of the function F(f*(x,t,s)), as e — > +0, where 
f*(x,t,s) S 7Y* is the approximation of the asymptotic distribution f*(x,t) £ £* . 

Then, for example, F(eS(x)) is defined as the weak asymptotics of the approximating function 

■fO. Taking into account the estimate for the function u and using the Lagrange 



F( W (f) 
theorem 



as e 



F(uj{z)) - F(0) = F'(Qw(z))u(z), < 6 < 1, 



we have 



\F(to(z))~F(0)\<K(l + \z\)- 3 . 

Thus, for the function F(iv(z)) — F(0) there exists an estimate analogous to that for the function ui(z) 
approximating the asymptotic distribution eS(x). Therefore, applying f(u>(^ )^ to a test function and 
performing the change of variables x — erj, we obtain 



J(e) = (F(0),ip(x)) 
= (F(0),<p.(x))+e 



F 



- F(0)1 <p(eri) dr) 



F 



(u(ri)) - F(0)] d?7 + 0(e 2 ), 



= (F(0),(^(x))+^(0) 
Thus, in the weak sense, we have 

F(eS(x)) = F(0) + sA5(x) + O v >{e 2 ) 

where the constant A = J F(u{rfj) — F(0) drj. 



+0, 



(2.29) 
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3 Infinitely narrow solitons to the Hopf and KdV 
equations 

1. The Hopf equation: u t + (u 2 ) x — 0. To find an asymptotic 0-p'(£ 2 ) infinitely narrow J-soliton-type 
solution to the Hopf equation ( |l.4| ) of the form 



i*(x,t) = uo(x,t) + g(t)e5(x - <f>(t)), e — ► +0, 



(3.30) 



where u (x,t) G C°°(R ), g(t),</>{t) G C°°(R) are the desired functions, one needs to substitute into the 



equation its approximation from 7i* of the form: 

u*(x, t, e) — uq(x, t) + g(t)eS(x — 4>(t), e), e > 0, 



(3.31) 



where, according to (2.25), the function eS(x, e) = w(f ) is used as an approximation of the asymptotic 
distribution eS(x). 

Using the first formula ( 2.28j ) for the squared asymptotic distribution ed(x), we obtain 



[u* £ (x,t)} 2 = ul + 2u gsS(x - <p(t)) + eg 2 flS{x - <p(t)) + V '(s 2 
Then, substituting U*(x, i) into the Hopf equation we have 

u *t + UK) 2 }* = u ot + (uq) x +sS(x - (/)) (^g t + 2gu 0x 
+ e5'(x - 4>{t)) ( - g<t> t + n g 2 + 2u .g) + O v , (e 2 ), e - 



-0. 



+o, 



where fl = j uj 2 (t]) drj. 

Using the well-known equality 



a{x)5\x) = a(0)8'(x) - a'(0)6(x), 



(3.32) 



(3.33) 



and equating the coefficients of e , ed and eS' with zero, we obtain the necessary and sufficient conditions 
for the right-hand side of equation ( |3.32| ) to be of the order of Ov{s 2 )'- 



u ot + (ul) x = 0, 



4>t - flg(t) - 2u (x,t) 



= 0, 

X=0(t) 

g t (t) = o. 



(3.34) 



In the case of constant background uq = const we have <j){t) = vt + 4>q, where v — Qg(0) + 2uq and (f>o 
is a constant which has the meaning of the coordinate of the initial position of the soliton. 

It is clear that the condition g t = 0, when uq =^ const, contradicts the physical intuition and the 
well-kno wn r esults about soliton behaviour. However, formally, such a structure can also exist, since 
system ( 3.3S ) is well-defined, but we can try to repair the situation by adding a new term which, after 
differentiation, has th e fo rm e5(x — ip). Therefore, we shall seek an asymptotic solution to t he Hopf 
equation of the form (L8) by substituting into the equation the smooth ansatz of the type (1.9): 



u*(x,t,e) = u (x,t) + g{t)uj(a(t)- — — ) + e(x, t)euj {- — — ; 



>0, 



(3.35) 



where a{i) G C°°, f2i = J uj{rf)dri, lu(z) is a S-type function which has a compact support or decreases 
sufficiently rapidly as \z\ — > oo, wo( j ) is an approximation of the Heaviside function. 
The weak asymptotics ( p. 35 ) has the form 



u*(x, t) = u (x, t) + g(t)n 1 eS(a(t)(x - <f)(t))) + e(x, t)e8(x - <j>(t)), 



-0. 



(3.36) 



This representation coincide with one that can be obtained when considering the smooth asymptotics of 
the asymptotic solution to the KdV equation with a small dispersion ~ e 2 [[j] , j| , in the weak sense as 
e -> +0. 
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In this case, according to Section ^, 

[u* e (x, t)] 2 = u 2 (x, t) + ^jj (2«o(a?, t)g(t)Qi + g 2 (t)n 2 ^j eS(x - <£(*)) 

+ 2u (x,t)e(x,t)e6(x - <p(t)) + O v >{e 2 ) 1 e -> +0, (3.37) 

where fii = J oj{r f) dg, 0,2 — J u> 2 (g) dg. 

Substituting ( |3.35 ) into the Hopf equation and using (3.3E), (3.37), up to the terms of the order of 
Oc(e 2 ), we obtain a relation analogous to ( |3.32 ). Then, setting the coefficients of e°, e0(x — </>(£)), 
eS(x — (j>(t)) and ed'(x — <p(t)) equal to zero, we arrive at the following result. 



Theorem 3.1 Suppose that for t € [0, T] there exists a smooth solution Uo(x,t) to the Hopf equation 

with the smooth initial condition uo(x,t) = Uq(x). Then, on the closed interval [0, T], the Hopf 

t—o . . 

equation has a solution in the form of an infinitely narrow 6-soliton (3.36) if and only if the unknown 

smooth functions Uo(x,t), g(t), e(x,t), a{t), <p(t) satisfy the following system of equations: 

u ot + (uq) x = 0, 
0t-2uo(0(t),t)-^(*) = 0, 



e(<f>(t),t) 



n? (g(t)/a(t)) t 



n 2 g{t) 
(e t (x, t) + 2(u (x, t)e(x, t)) x ) 



X>(j>(t) 



= 0, 
= 0, 



(3.38) 



where Q,\ = Jaj(n)dri, VI2 — j oj 2 (n) dn . 



Remark. System ( |1.12|) s tudied in the Introduction can be obtained from system ( 3.38 ) if we take 
the solution of equation (1.20) as the approximation uj{a{t)(x — ip(t))/e). In this case we have 



uj = cosh 2 (t), 



y/g(t)/6 



and system ( ^.3§| ) turns into systems ( |l.l2| ), where the number of unknown functions is equal to the 
number of equations. 

2. The KdV equation: u t + (u 2 ) x +e 2 u xxx — 0. This equation has the exact one-soliton solution (|l.2|) 



^ 2 (\/f (l -i 9f)/e 



u(x,t,e) = g 



whe re g = const is the amplitude of the soliton |24| . The weak asymptotics of this solution has the form 

O) ' 

2 

u(x, t) = geS(x — -^gt), e ~ * +0- 

Let us study the dynamics of propagation of an infinitely narrow deformed soliton solution of the 
KdV equation. To this end, we consider a smooth ansatz of the form 

x — (h(t) x — <b(i) 
u*(x,t,E) = u (x,t)+g(t)w(a(t) — ) + e(x, t)eu>n{ — ), e > 0, (3.39) 

£ E 



where a(t) = sj ^i, u{g) = ch z (g), uj (|) = 1 — w o(f )■ ancl , according to (2.26), 



X 



LUi(g)dg. 



Here Wo(f) an d oj q (|) are approximations of the Heaviside functions 6{x) and 9(—x), respectively, 
Juj 1 (g)dg = l. 

The weak asymptotics of the right-hand side of ( |3.3S| ) has the form ( 3.36 ). 
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According to Definition 1.1, the infinitely narrow soliton-type solution of the KdV equation is defined 
as the solution of the following system of the two equations 



L KdV [u] =u t + (u 2 ) x + e 2 u xxx = Ov{e 2 ) 
uL Kdv [u] = (u 2 ) t + ^{u 3 ) x + e 2 2uu xxx = Ov{s 2 ). I ' > ' 4U1 

Let us write the dispersion term of the second equation in the form Q[u, e) = e 2 2uu xxx — e 2 [{u 2 ) xx — 
i{u x ) 2 ] x . 

Up to the terms whose asymptotics are of the order of Ox>'(e 2 ) in V we obtain 

[u*(x,t,e)} 2 = u 2 , + g 2 iv 2 + 2u gu + 2u eeuj- , 
[u*(x,t,e)} 3 = + g 3 LU 3 + 3uog 2 uj 2 + Sulguj + 3uleeuj^ . 



Then for all <p(x) S V: 



u(a(t f ^ )y(g) dx = ^t)J u (v)<p(<f>(t) + -^v) dn = ^y^WW) + °( £2 )' 



where Q± — fu>(r))dr) — Jch. 2 (r))drj = 2. Therefore, 

Ox 



V e J ait) 

In a similar way we find that 



UJ 



^ a{t) x _lM^^ e5{x . m ) + OMe 2 l 



where tt 2 = J ^ 2 (v) dr) = J ch" 4 (r?) <2?7 = §, fl 3 = J u 3 {v) dr] = J ch" 6 ^) cfy = if. 

Now we find the asymptotics of the dispersion term Q(u*(x,t,e),e). It immediately follows from 



Q3.37D that 

g(u*(x, t, e),e) = e 2 [{{u*) 2 ) xx - 3(u* x ) 2 ] x = -ie 2 [{ul) 2 ] x + O v >{e 2 ). 



^From (3.39) we have 

u%{x, t, e) = u 0x + guj x + e x sujg + eeuj^ x . 

Squaring this expression and analyzing its components, we see that up to the terms of the order of 
Ox>'(e 2 ) we have 

g(u*{x,t,s),e) = -3e 2 g 2 [(Lu x ) 2 ] x + Ov(s 2 ). 
After calculating the weak asymptotics of this expression, we obtain for all f(x) G T> 



{g(u*(x,t,e),e),<p(S)) =-3g 2 (t)e< 



<p{x) dx 

= -3g 2 (t)ea(t) I [u/( V )] =VW(t) + — V ) dr) = -3n 4 g 2 (t)a(t)e^(t)) + 0(e 2 ), 
J ot{t) 

where il 4 = JV(?7)] 2 dr] = 4 / ch" 6 (r?) sh 2 (r?) dr] = if. 

It follows that the weak asymptotics of the dispersion term is 

G(u* e (x,t,),e) = -3g 2 (t)a{t)n 4 e5'(x - c/>(t)) + V '(e 2 ). (3.42) 

Substituting the asymptotics for u, oj 2 , uj 3 obtained above into ( 3.39 )-( j3.4l| ) we have 



u*Jx, t, )=u + g — e5(x - <j>(t)) + e(x,t)e6(-x + 4>(t)) + Ov>{e z ), 
a 

:(x,i)] 2 = ul + {g 2 — + 2u g — )eS(x - <p(t)) + 2u ee8{~x + <j>(t)) + O v ,(e 2 ), 
la a j 
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[u*Jx, t)] s =u 3 + (g 3 — + 3u g 2 — + 3u 2 g— )eS(x - <£(*)) + 3u 2 ee6(-x + 0(t)) + 0^(e 2 ). 

After the substitution of the coefficients ilk, k = 1,2, 3, into these expressions we obtain 

u*(x, i) = w (a;, t) + 2V6g 1/2 {t)eS(x - 0(f)) 

+e(x,t)s8(-x + <j>(t)) +Ov>(e 2 ), 

[u* e (x,t, )] 2 = u 2 (x,t) + V6{^g 3/2 (t) + 4«o(0(t), t)ff 1/2 (*)}^(-a: + 0(t)) 

+2u o (a?,t)e(a: ) i)eff(-a; + 0(t)) + 2 >#(e 2 ), (3.43) 

[«*(*, i)] 3 = u 3 (x, t) + V&{ ^.9 5/2 (i) + 4uo(0(t), f)<? 3/2 (i) 

+6^(^(t),*)s 1/a («)W(t),*)}^(a; - 0(t)) 

+3w 2 ,(a;, t)e(ar, <)e6»(-a; + </>(*)) + Ox>/(£ 2 ). 

According to Definition [Tl], substituting asymptotics ( 3.43 ) into the first equation of system ( |3.40| ) 
and setting the coefficients of e , £(5, £<$' and sO equal to zero, we obtain the necessary and sufficient 
conditions for the right-hand side of the equation LkcIv[u} — to be of order Ov(s 2 )- 



,0 2 



f2i 



eS' 
eS 

e6 



2 , 



= o, 

-<f>t + 2u (4>{t),t) + lg{t) = 0, 



(3.44) 



(e t (x, t) + 2(u (x, t)e(x, t)) x ) 



x<<j>(t) 



The second and third equations of system ( |3.44| ) imply the equation for the jump of the amplitude of 
a small shock wave: 



e(<p(t),t) 



2 gV 2 (t) 



(3.45) 



Following Definition 



1.1 



we now substitute asymptotics (3.42) into the second equation of system 
(3.4C) and, setting the coefficients of e°, eS, eS' and e9 equal to zero, we obtain the necessary and 
sufficient conditions for the right-hand side of the equation uLkov [u] = to be of order Ox>'(e 2 ): 



6° 

e5' 
eS 



e6 : 



K)t + K w o)x 



0, 



■0t+2«oWt),t) + |ff(t))(fl(t) + 3uo(^(t),t)) = 0, 
V6g(t)(g(t) +2u ((t>(t),tf 



-§«o(0(t), <).9 3/2 W ( - e{^{t), t) - |V6 
((2uo(a, t)e(x, t)) t + i(u 2 (x, t)e(x, t)). 



gt{t) 
g 3/2 W 



x<<j>(t) 



(3.46) 



= 0, 
= 0. 



Since uo(x,t) and e(x, t) are smooth functions, the first and the last equations from systems Q3.44 ) 
and ( |3.46[ ) are equivalent. It follows from the third equation of system (3.46) and equation ( 3.45| ) that 
g(t) + 2u ((j)(t ). t) = co nst. 

Relations ( 3.44 ) - ( 3.46 ) imply the following theorem which determines the dynamics of a single 
deformed soliton to the KdV equation. 

Theorem 3.2 Let us assume that for t € [0, T] there exists a smooth solution Uo(x,t) to the Hopf 

equation with the smooth initial condition uo(x,t) = u9.(x). 

t—o 

Then the KdV equation on the closed interval [0, T], up to Ov{£ 2 ), has an infinitely narrow S-soliton- 
type solution 

u*(x,t) = Uo(x,t) + 2V6g 1/2 {t)eS(x - 4>{t)) + e(x,t)e0(-x + <j>(t)), 
if and only if the unknown smooth functions Uo(x,t), g(t), e{x,t), <j)(t) satisfy the system of equations 



U 0t + (Uo) a 



= o, 



4> t = 2u (<f>(t),t) + h(t), 



(e t (x, t) + 2(u (x, t)e(x, t)) x 



x<4>(t) 



$(O) + 2u8(0(O)), 
0, 

3-/6 sM 



(3.47) 
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Thus system of equations coinsides with system (1.15) studied in the Introduction. 



4 Infinitely narrow solitons to the equations u t + 

(f(u)) x = and u t + (f(u)) x + £ 2 u xxx = 

1. The equation u t + (f(u)) x — 0. Consider the infinitely narrow soliton-type solution of this equation, 
substituting smooth ansatz (1.22) into the equation. This ansatz is a generalization of the smooth one 
( 3.3 5| ) and has, up to 0-d'(e 2 ), a weak asymptotics of the form (1.23). In this case, in (1.22) wo(§) is an 
approximation of the Heaviside function 6{x). 

By analogy to ( [2.29 ) , define the asymptotic distribution 



f(u (x, t) + Clx(t)eS(x) + e(x, t)e6(xfj 
as the weak asymptotics of the approximating function as s — * +0 

f(u (x, t) + e(x, t)eu (^)J . 

It is clear that 

f(u + u(^,t) +e£UJ {^)) = f(u + w(pt)) + e/'(uo)ewo(^) + 0(e 2 ). 
In addition, according to the Lagrange theorem, we have 

/(n + LJ^t)) - /(n ) = /' (uo + 0u,(f 



where < G < 1. Hence, for the function f\UQ + u>(j,t)j — f(uo) we have the same estimation with 
respect to r as for the function ui(t, t), 

Applying the function f(uo + cj(|,i) + eewo(f )J to a test function and performing the change of 
variables x = et, we obtain, as in deduction of formula ( 2.29 ), 



(X X \ 

u (x, t) + u(-,t) + esujo(-) 



e e ■ 

= f(u (x,t)) + A(0,t)e5(x) + f(u o (x,t))e(x,t)s0(x) + 0^(e 2 ), 
where the function 



A(x,t)= / f(u (x,t)+u(T,t))-f(u (x,t)) 



dr. 



(4.48) 



(4.49) 



Substituting u*(x, t), in the form of singular ansatz ( 1.23| ), into the initial equation and using ( |4.48| ) 
we obtain 

(/«)), = u ot + (f(u )) x + { - &(t)fii(t) + A(c)>(t),t)}e6'(x - </>(t)) 



+{(Oi(*))t - e(x, t)(f> t (t) + e(x, t)f'(u (x, t))}eS(x - (f>(t)) 
+{e t (x,t) + (f'(u a (x,t))e(x,t)) x }e9(x - </>(t)) + O v >(e 2 ), e - 



-0. 



Using formula ( |3.33] ) and setting the coefficients of e°, e8{x — 4>), e5'(x — </>), e9{x 
zero, we arrive at the following result. 



equal to 



Theorem 4.1 Suppose that for t € [0, T] there exists a smooth solution uo(x,t) to equation uot + 

(f( u o))x — with the smooth initial condition Uo(x,t) = u9,(x), where f(u) is a smooth function. 

t=o 
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Then, for t G [0, T], the initial equation has the infinitely narrow S-soliton solution (1.22) if and only 
if the unknown smooth functions Uo(x,t), e(x,t), 4>{t) satisfy the following system of equations: 



(«i) t + f'(u (<f>(t),t))- 



v -™ ni(t) 

e t (x,t) + (f'(u (x,t))e(x,t)) 



uot + (f(uo)) x = 0, 

0tni(t) + A(^(i),t) = o, 

= o, 

x>tp(t) 



(4.50) 



where £li(t) = J u>(r,i) dr and the function A(x,t) is defined in (4.49). 

It is easy to verify that if the inequality f"(u) > is satisfied, then for fii(i) > system (4.54) has 
exactly the same properties as system (1.14). 
2. The equation u t + (f(u)) x + £ 2 u xxx — 0. 

As it was said in Introduction, for some convex smooth f(u) this equation has exact one-soliton 
solutions. In this case we use our approach to describe the dynamics of propagation of a deformed 
infinitely narrow soliton. By virtue of Definition 1.1, as in the case of the KdV equation, the infinitely 
narrow soliton- type solution (1.23) is defined as a solution of the two equations 



L[u] =u t + {f{u)) x + e 2 u xxx 



Or 



uL[u] = {u 2 ) t + (f(u)) x + e 2 2uu xxx = V '{e 2 ), 



(4.51) 



where f(u) = 2uf{u) — 2 J™ f(z)dz. Note that the ansatz ( 1.23| ) has the approximation (1.22), where 
ujq = luq (-) is the approximation of the Heaviside function 6(—x). 

The dispersion term in the second conservation law is just the same as for the KdV equation: 

G(u,e) = e 2 2uu xxx = e 2 [(u 2 ) xx - 3(u x ) 2 ] x , 



and its weak asymptotics has the for m (3.42 ) and fl^i) = J[ui' t (t, t)] 2 dr. 

Substituting the smooth ansatz ( 1.22j ) into the initial equation, just as for the KdV equation from 
Section ||.2, we obtain, 

[u* e {x, t)} 2 =u 2 + |n a (t) + 2u n 1 {t)}e5(x - <j>(t)) + 2u ee8(~x + 4>{t)) + V '(s 2 ), 

where Qi( t) = J u{T, t)dT Vl 2 (t) = J lu 2 (t t)dr . 

From jL22] ) and (jL2|), using formula ([4.48[), we find 



f(u (x,t) + Cli(t)eS(x) + e(x,t)e6(-x)) 

= f(u (x,t)) + A(fl,t)eS(x) + f'(u (x,t))e(x,t)e9(-x) + Ojy(e), 
where f'(u) — 2uf'(u) and 



-0, 



A(x,t)= / f(u (x,t)+uj(T,t))-f(uo{x,t)) 



dr. 



(4.52) 



(4.53) 



Obviously, the system of equations derived from the first equation in (4.51) coincides with system 
( 4. 50] ) from Theorem 4.1 which describes the dynamics of the formal soliton solution of the equation 
ut + {f{u)) x = 0. 

We find the system of equations which follows from the second conservation law uL[u] = in (4.51). 
To this end, substitute the asymptotics [u*(x,t)] 2 , ( 4.52| ) and (3.42) into the second equation (4.51) and 
set the coefficients of e°, eS, eS' and e9 equal to zero. 

Setting the coefficients of e° and e9 in the obtained system equal to zero, we have the equations 



(«2)t + (/M) 

(2u e) t + (2u f'(u )e 



x<4>(t) 



which, in virtue of smoothness of the functions uo(x,t), e(x,t), f(u), coincide with the corresponding 
equations of system (4.50) of the first conservative law L[u] — 0. 
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After setting the coefficient of e5 equal to zero, we find 

(n 2 (t) + 2«o(^(t), (*)) t + 2u {4>{t), t)e(<f>(t),t) (/'K) - 0t(t) 



0. 



Using the second and the third equations of system (4.50), we bring the last equation to the form 

(4.54) 



j t n 2 (t) +2Q 1 (t)j t u ( ( f>(t),t) =0, 



which coincides with the third equation from (3.47) for the KdV equation when a(t) = y 2^1, 
Consider the coefficient of eS' and prove 



Lemma 4.1 The coefficients of sS' in the conservation laws L[u] = and uL[u] = (1-51) coincide. 



Proof. Setting the coefficient of eS' in the left-hand side of the equation uL[u] = equal to zero, we 



find: 



0t(t)(n 2 (*) + 2uofa(t),t)fii(t)) +M4>(t),t) - 3n A (t) = o. 



Let m(x, i) = uq(x, f) + u>(t, f) be an exact soliton solution to the equation ut + {f(u)) x + e 2 u x 



where r 



Substitute u(x, f) into this equation and consider the coefficient of p 
- <f> t w r + (f(u + (jj)) T + U) TTT = 0. 



(4.55) 

•:x — 0, 



(4.56) 



Integrating equation ( 4.56| ) with respect to r from — oo to r and taking into account that ui(r,t) and 
its derivatives with respect to r tend to zero as |r| — > oo, we obtain 



f W + f(u + u) + LO T T = f(u ). 



(4.57) 



Here the integration constant f(uo) has been derived from the boundary conditions. 

Multiplying equation ( 4.57[ ) by ui T and integrating it again with respect to r from — oo to r, we have 



4>tU 2 + 2 / f(uo + w)w T dr — 2/(uq) / w r dr + (w T ) 2 = 0. 



(4.58) 



Let us introduce the function /i(it) = L f(z) dz, /i(0) = 0. 
Since ^ 



and since cj(t, i) — * for r — > — oo, 

/(ito + w)w T dr 



d - 



dr 



/i (u + uj) dr = /i (it + w) - /i (u ) 



Now equation (4.5E) can be rewritten in the form 



~(j) t io 2 + 2 [/x («o + u) - fx (uo)J - 2/(tto)w + K) 2 = 0. 
Integrating this equation with respect to r from — oo to oo, we find 

- 4> t n 2 - 2f(u )n 1 + 2Ai + fi A = 0, 
where Qi(i) = / u(r,t)dr, Sl 2 (t) = J w 2 (r,i)dr, fi A (f) = / (uJ T ) 2 (T,t) dr, 



Ai(x,t)= / /i(u (x,t) +w(r,t)) - /i(uo(i,t)) 



dr. 



(4.59) 



(4.60) 



Multiplying (4.57) by cj and integrating the obtained equation with respect to r from — oo to oo, we 



find 



6(^2 + 



f(u +uj) - f(u ) 



lu dr — f2 A = 0, 



(4.61) 
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where the last term has been derived by integrating the expression J co tt lu dr by parts. 
Subtracting equation ( ft.59|) from the doubled equation ( 4.61 ), we have 



(p t O 2 + 2f(u )O 1 +2 f(u +Lj)-f(u ) 



udr- 2Ai - 3fi A = 0. 



(4.62) 



Since /(it) = 2uf(u) — 2/i(it), where fi(u) = J™ f(z)dz, the function A(x,t) from (4.53) can be 
rewritten in the form 



A(x,t) 



(ua + uj)f(u +uj)- u f(u ) 



dr -2 



f(u Q +uj) - f(u ) 



dr. 



Adding and subtracting /(wq)w in the first integral and comparing the obtained expression to formulae 



(4.48) and (4.6C) for A(cc,t) and Ai(x,t), respectively, we obtain 

A(x, t) = 2w A(x, t) + 2/(u )fti + 2 \f(u + w) - f{u ) uodr- 2Ai(x, t). 



(4.63) 



Setting the coefficient of eS' in the first conservation law L[u] = equal to zero and taking into 
account Theorem 4.1, we arrive at the equation 



-&fii(t)+A(#f),t) =0. 



Substituting this relation and relation ( 4.63 ) into equation ( 4.62 ), we obtain equation ( 4.55| ), that is, 
the coefficient of sS' in the second conservation law = 0. Thus, the equations obtained by setting 

the coefficients of ed' in both conservation laws ( 4.51 ) equal to zero are equivalent. 

The proof of the lemma is complete. 

As in Section ||.2, this implies the theorem which determines the dynamics of propagation of the 
soliton solution of equation under consideration. 

Theorem 4.2 Suppose that for t € [0, T] there exists a smooth solution uo(x,t) to the equation Ut + 
(f{u)) x = with the smooth initial condition uo(x,t) = Uq(x). Suppose that the problem (1.21) has a 

i=0 

solution and f"(u) > for u € [wi,^]- 

Then the e quation ut + (f(u)) x + s 2 u xxx = for t € [0, T] has an infinitely narrow delta- soliton-type 
solution ( 1.23 ) if and only if the unknown smooth functions uo(x,t), g(t), e(x,t), <f>(t) satisfy the 
following system of equations: 



uot + (f{uo)) x = 

-&fii(t)+A(#t),t) = 

jn 2 (t) + 2n 1 {t)ju m),t) = o 
(nO.+ f/'MW),*)) k{mt] 



(4.64) 



e t (x,t) + (f'(u (x,t))e(x,t)) 



e(4>(t),t) = 0, 
= 0, 



x<<p(t) 



where ill = J u>(r,t) dr, 0,2 — Jw 2 (t, t)dr, 

A(<f>(t),t) = [ [f(u Q ( ( f>(t),t)+g(t)co(T,t)) - /(uo(0(t),t))] dr. 



It is easy to see that the system given in Theorem 4.2 and system (1.15) discussed in the Introduction 
are quite similar. 

In conclusion, we note that all constructions carred out in this paper are formal and the weak asymp- 
totics method is still not completely developed. 

This work was partially supported by the Russian Foundation for Basic Research (Grant Nos. 97-01- 
01123 and 99-01-01074). 
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